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ABSTRACT 

A  new  approach  is  described  to  the  connection  of  wave  amplitudes  across 
the  turning  points  and  singular  points  of  second-order,  linear,  analytic, 
ordinary  differential  equations  which  can  describe  the  modulation  of  physical 
waves  or  oscillators.  The  general  class  of  singular  points  thereby  defined 
(Section  2)  contains  many  irregular  ones  of  greeter  complexity  than  have  been 
accessible  before;  however,  genuine  coalescence  of  singular  points  is  not  here 
considered.  The  asymptotic  connection  formulae  are  shown  to  result  directly 
from  the  branch  structure  of  the  singular  point  ('action  3),*  indeed,  to  a 
first  approximation,  they  reflect  merely  the  gross,  local  branch  structure. 

The  proof  (Section  4)/  relates  the  local  structure  of  the  solutions  at  the 
singular  point  to  the  asymptotic  wave  structure  by  a  limit  process  justified 
by  symmetry  bounds . 
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SIGNIFICANCE  AND  EXPLANATION 


This  work  concerns  the  modulation  of  waves  or  oscillating  systems,  which 
pervade  all  the  science  and  engineering  disciplines.  Modulation  occurs  when 
waves  travel  through  an  inhomogeneous  material  in  which  the  local  propagation 
velocity  differs  from  place  to  place,  but  the  differences  are  small  over  a 
distance  of  only  a  wavelength  —  a  very  common  case  in  the  sciences  and 
engineering.  The  resulting  change  to  the  waves  is  mostly  gradual,  but 
occasionally  drastic,  as  at  a  shadow-boundary,  where  oscillation  turns  into 
decay  and  quiescence  over  just  a  few  wavelengths.  When  this  phenomenon  can  be 
analyzed  via  an  ordinary  differential  equation,  such  a  boundary  is  called  a 
transition  point. 

At  first,  only  the  simplest  transition  points  representing  the  most 
typical  shadow  boundaries  were  studied.  But  then  some  phenomena,  such  a  wave 
reflection  and  scattering  cross-sections,  came  to  be  traced  to  hidden 
transition  points  that  become  visible  only  when  real  distance  (or  time)  is 
embedded  in  its  complex  plane.  When  the  material  properties  vary  in  a  general 
manner,  (which  can  often  be  observed  only  incompletely)  the  hidden  *-ransition 
points  can  have  arbitrarily  complex  structure.  The  following  v_  agents  a 
new,  more  direct  and  more  general  approach  to  the  connection  of  waves  and 
shadows  across  transition  boundaries.  It  aims  to  furnish  a  basis  for  more 


CONNECTION  FOR  WAVE  MODULATION 


R.  E.  Meyer  and  J.  F.  Painter 
1.  Introduction 
The  semi-classical  Schroedinger  equation 
(1)  e^d^w/dz^  +  p  w(z)  *  0 

with  small  parameter  e  and  analytic  coefficient  function  p(z)  is  central 
to  a  vast  class  of  oscillation  and  wave  modulation  problems  in  physics  and 
other  sciences.  Particular  interest,  especially  for  scattering  theory, 
attaches  to  the  "WKB"  problem  of  connecting  the  wave-approximations  to 
solutions  across  roots  or  singular  points  of  p(z).  The  following  introduces 
a  connection  method  which  is  simpler  and  more  general  than  any  advanced  before 
[Zwaan  1919,  Langer,  1931,  Painter  and  Meyer  1981].  Simplification  and 
clarification  of  connection  theory  is,  in  fact,  the  whole  objective  of  the 
study  to  be  reported,  and  generalization  was  used  only  as  a  help  towards  it. 

One  reason  why  this  objective  has  proved  elusive  over  the  generations  may 
be  that  the  general,  second-order,  linear  differential  equation,  of  which  (1) 
is  the  normal  form,  encompasses  too  many  disparate  phenomena.  The  present 
study  focuses  on  only  those  forms  of  (1)  which  are  genuine  Schroedinger 
equations  in  the  sense  that  they  can  describe  the  modulation  of  physical  waves 
or  oscillators.  This  subclass  is  characterized  in  Section  2  in  terms  of  its 
admissible  (turning-point  and)  singular-point  structure.  To  attempt  only  one 
step  at  a  time,  moreover,  genuine  coalescence  of  singular  points  is  excluded. 
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This  leaves  a  large  class  of  singular  points,  none  the  less,  because  the 
potential  functions  p(z)  of  (1)  in  the  sciences  must  be  defined,  it  not  by 
speculation,  then  by  measurement,  in  which  case  they  can  be  known  only 
imperfectly.  In  addition,  it  has  long  been  recognized  that  scattering 
matrices  may  depend  greatly  on  singular  points  of  p(z)  away  from  the  real 
axis  of  time  or  space,  in  which  case  there  are  scant  physical  grounds  for 
restricting  their  nastiness.  The  characterization  of  p(z)  cannot  therefore 
be  very  specific  and  must  include  arbitrarily  irregular  points  of  (1) . 
Certainly,  multi-valued  functions  p(z)  must  be  the  norm,  rather  than  the 
exception.  All  the  same,  modulation  implies  a  certain  structure  (Section  2). 

The  multi-valuedness  of  p(z)  in  (1)  implies  that  the  solutions  w(z) 
must  normally  be  multi-valued,  and  the  main  thesis  to  be  propounded  is  that 
this  multi-valuedness  is  the  source  of  the  connection  problem  and  that  the 
asymptotic  connection  formulae  solving  it  are  a  direct  manifestation  of  the 
branch  structure  of  the  potential  p(z)  at  the  singular  point.  This 
fundamental  view  of  connection  is  adopted  also  by  Olver  [1974,  pp.  481,  482] 
for  regular  and  isolated  singular  points.  Our  objective  is  to  show  how  it  can 
be  extended  to  large  classes  of  very  irregular  ones  and  what  new  insights  into 
the  nature  of  connection  emerge  therefrom. 

To  obtain  ammunition  for  this  thesis,  local  solution  representations  near 
the  singular  point  have  been  developed  in  a  companion  paper  [Meyer  and  Painter 
1982,  hereafter  referred  to  as  IPM]  and  are  summarized  in  Section  3.  They 
focus  on  a  particular  fundamental  system  yg ,  ym  of  (1)  in  which  ym  has  a 
milder  singularity  than  yg  which,  in  turn,  contains  no  additive  multiple  of 
ym .  This  makes  the  pair  a  characteristic  representation  of  the  branch 
structure  of  the  singular  point.  The  representation  is  a  local  one,  in  the 
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first  place,  but  turns  out  to  have  a  striking  two-scale  structure:  in  the 
framework  of  a  natural,  independent  variable  x  (Section  2),  the  solutions  do 
not  depend  on  the  complex  parameter  e  in  (1)  separately,  but  only  on  the 
variables  x  and  ex.  A  key  property  of  the  representations  is  that  they  are 
global  in  the  oscillation  variable  x,  even  though  only  local  in  the 
modulation  variable  ex,  and  can  be  extended  to  similar  bounds  on  a  loss  of 
symmetry  of  the  solutions  with  increasing  distance  from  the  singular  point. 

To  document  the  thesis,  these  symmetry  bounds  are  shown  in  Section  4  to 
adnit  a  class  of  limit  processes  in  which  both  |x|  +  *  and  |ex|  ♦  0  and 
furthermore,  the  symmetry  of  yg  and  ym  is  essentially  preserved.  In  other 
words,  asymptotic  approximation  of  WKB-type  characterized  by  dominance  and 
recessivity  is  shown  to  become  available  before  local  structure  has  been 
lost.  Existence  of  such  limit  processes  translates  immediately  local 
information  on  the  structure  of  yg  and  ym  at  the  singular  point  into 
information  on  the  multivaluedness  of  their  asymptotic  wave-represenation . 

But,  the  latter  information  is  what  connection  theory  seeks. 

To  make  our  point,  it  will  be  sufficient  to  document  it  in  the  generic 
case,  leaving  aside  the  "Frobenius  exceptions"  [Olver  1974  pp.  150]  which 
involve  logarithmic  branch  points  for  regular  points  and  consequent  loss  of 
the  symmetry  bounds,  in  the  general  case.  The  validity  of  the  final 
connection  formulae  (Section  4)  for  the  exceptional  cases  is  strongly 
indicated  by  the  results  of  Painter  and  Meyer  [1981],  but  the  approach  adopted 
in  the  following  seems  unsuited  for  a  simple  proof.  It  is  also  sufficient  to 
document  our  point  at  the  instance  of  the  first  asymptotic  approximation , 
which  connects  the  wave-amplitudes,  and  attention  will  be  restricted  to  it.  A 
more  exhaustive  description  may  be  thought  desirable  for  the  sake  of 
completeness,  and  most  of  all,  error  bounds  are  desirable.  It  may  be  noted 


that  the  representations  used  are  obtained  by  the  standard  method  of  Volterra 
integral  equations,  which  is  precisely  the  method  leading  to  effective  error 
hounds  [Olver  1974].  The  sketch  of  the  representation  method  given  in  Section 
3  makes  clear  that  such  bounds  would  emerge  in  terras  of  pointwise  and  integral 
bounds  on  a  certain  irregularity  function  $(£x)  arising  in  the 
characterization  of  the  singular  point  (section  2).  In  the  general  case, 
however,  that  function  is  barely  specified  and  can  tend  to  zero  with  ex 
arbitrarily  slowly,  so  that  it  appears  doubtful  whether  the  bounds  would  give 
much  satisfaction  to  the  numerical  analyst.  Most  of  all,  however,  we  suspect 
that  the  present  proof  of  connection  may  not  remain  the  simplest  one  for  long, 
in  which  case  present  attention  to  error  bounds  and  higher  approximations  may 
be  premature . 


2.  Modulation  Equations 


Equation  ( 1 ) 

2  2  2 

£  d  w/dz  +  p( z ; £ )w( z )  -  0 

is  one  of  a  large  family  of  normal  forms  of  the  general,  linear,  second-order 
differential  equation  and  constructive  general  statements  are  difficult  in 
such  an  Indefinite  frame .  By  oonstrast, 

dx/dz  ■  ip^2  /£ 

defines  the  Lionvilie-Green  or  VKB  or  Langer  variables  x  and  ex  based  on 
the  local  wavelength  or  period,  which  have  long  been  recognized  as  the  natural 
ones  for  wave  modulation.  Physical  specifications,  e.g.  for  scattering, 
relate  directly  to  them  and  if  z  differs  substantially  from  x,  it  can  at 
best  measure  distance  in  legal  units .  The  natural  formulation  of  physical 
problems  of  wave  modulation  is  therefore  in  terms  of  x  or  £x,  from  the 


start,  which  will  avoid  the  extraneous  difficulty  of  describing  the  global 


transformation  between  ex  and  z,  which  has  no  physical  significance  and 


can  be  a  very  complicated,  multi-valued  map. 

The  exclusion  of  coalescence,  in  order  to  confine  attention  to  one 
singular  point  at  a  time,  restricts  not  their  total  number,  but  only  how  fast 
they  can  approach  each  other  as  e  -*■  0 .  When  this  is  not  fast  enough  to 
introduce  genuine  coalescence,  a  rescaling  [Meyer  and  Guay  1974]  permits  the 
elimination  of  the  main  e-dependence  from  lt(z;£)  and  therefore,  not  much 
generality  is  lost  by  ignoring  the  residual  e-dependence,  for  simplicity, 
p(z)  will  accordingly  be  taken  independent  of  e  in  what  follows. 

The  main  property  distinguishing  the  wave  modulation  equations  among  the 
larger  class  of  equations  (1)  is  that  the  natural  variable  x  must  be 
definable ,  for  otherwise,  not  even  the  concepts  of  wavelength  or  period  could 

V, 

exist  for  (1).  An  additional  requirement  arises  as  follows.  If  p  ^  be  non- 
Integrable  at  a  singular  point,  then  that  point  is  seen  to  correspond  to  no 
x  6  (C  and  hence,  represents  not  a  genuine  singularity  but  a  device  for  re¬ 
interpreting  radiation  conditions  as  a  singular  point  in  the  z-planes .  Such  a 
device  has  been  used  at  times  in  quantum  mechanics,  but  is  excluded  here  to 
concentrate  on  the  class  of  genuine  singularities  of  modulation.  For  that 
class,  the  singular  point  of  (1)  must  correspond  to  a  definite  point  x. 
Without  loss  of  generality,  both  will  be  identified  with  this  origin. 

For  an  effective  formulation  of  this  notion  of  the  most  general  wave 

modulation  equation  (short  of  coalescence),  it  should  be  expressed  in  terms  of 

the  natural  variable  x.  Accordingly,  the  following  is  based  on  the  premise 

1/. 

that  a  branch  r(x)  of_  p4  is  definable  as  an  analytic  function  on  a 
punctured  neighborhood  of  x  =  0  which  is  a  Riemann  surface  including  the 
interval  (-it ,2*0  of  arg  x  so  that 
(2)  i  dz/dx  =  e  r  ^ 
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is  inteqrable  at  x  =  0.  (In  conventional,  turning-point  terminology,  such  a 
Riemann  surface  element  comprises  three  adjacent  Stokes  sectors.) 

When  the  Schroedinger  equation  (1)  is  transformed  to  the  natural  variable 
by,  say,  simply  setting 

w(z)  =  y(x)  , 


it  takes  the  form 

(3)  y"  +  2r-1r *  y'  -  y  =  0  , 

which  shows  that  the  wave  development  is  controlled  by  the  "modulation 
function" 

-  V- 

r’(x)/r(x)  =  ( ie/2 ) d(p  '2  )/dz  , 

rather  than  by  the  potential  function  p(z)  directly.  This  illuminates  why 
it  has  long  been  known  that  the  singular  points  of  (1)  should  really  include 
the  roots  of  p(z)  ("turning  points").  It  is  also  seen  that  the  modulation 
function  has  a  particular  structure:  since  p(z)  =  r4  is  a  function  of  z 
independent  of  e,  it  follows  from  (2)  that  ex  is  also  such  a  function  and 
in  turn,  that  xr'/r  depends  on  x  and  e  only  through  the  product  ex.  A 
secondary  hypothesis  to  be  now  adopted,  because  it  simplifies  the  theory  of 
connection,  is  that  a  limt  of  xr’/r  as  ex  +  0  can  be  identified, 

(4)  xr'/r  Y  e  C  as  0 

uniformly  in  the  Riemann  surface  sector  A  of  ex  in  which  xr'/r  is 
defined  locally  near  ex  =  0.  These  two  hypotheses  also  define  the  framework 
of  the  analysis  of  [IPM] .  A  statement  equivalent  to  (4)  is  that  the  (fourth 
root  of  the)  potential  r(x)  can  be  written  in  the  form 

V,  y 

(5)  p  4  =  r ( x)  =  x  p(£x) 

where  P(S)  is  a  function  analytic  on  the  Riemann  surface  element  A  with 
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the  property 


(6) 


(£/P)dP/d£  =  :  $(£)  +  0  as  £  *  0 


uniformly  in  A,  because  <J> ( e x)  =  xr*/r  -  Y. 

To  make  the  structure  of  the  theory  more  readily  apparent,  it  will  help 
to  abbreviate  the  notation  by  the  convention  that  a  function  symbol  such  as 
g(x)  is  always  understood  to  denote  a  function  of  both  x  and  ex.  By 
constrast  a  Greek  symbol  such  as  'I'(C)  will  always  denote  a  function  of 
5  =  £x  only.  If  such  a  function  has  the  property  (6), 

(£/<!>) dij>/d£  >0  as  £  +  0 

uniformly  in  A,  then  it  will  be  called  mild;  it  implies  that  <|>(C) 
varies  near  £  =  0  less  than  any  non-zero  real  power  of  £: 

VV  >  0,  |£Ui|»±1(  +0  as  £  +  0  . 

In  particular,  the  limit  Y  postulated  in  (4)  is  thus  seen  from  (5)  to 
represent  the  exponent  of  the  "nearest  power"  of  x  in  the  (fourth  root  of 
the)  potential,  and  the  basic  integrability  premise  defining  physical 
Schroedinger  equations  implies 

Re  y  <  \  • 

The  general  class  of  singular  points  of  Schroedinger  equations  thus 
defined  includes  very  irregular  ones,  in  addition  to  all  the  turning  points 
covered  in  the  literature  (Painter  and  Meyer  1981).  For  banger's  [1931]  class 
of  fractional  turning  points, 

is  analytic  and  nonzero  at  z  =  0, 

*<£)  -  l  Y5M-2m 

n=1 

and  the  solutions  of  (1)  and  (3)  are  approximable  in  terms  of  Bessel  functions 
[Langer  1931,  Olver  1977].  For  othe<-  singular  points,  however,  no  simple 
approximands  in  terms  of  '‘'liw  . motions  appear  likely.  local 


approximations  have  been  constructed  in  [I PM]  to  provide  support  for  the 
present  study  even  under  the  vague  assumptions  just  sketched,  which  acknit 
functions  ♦  (£)  of  arbitrary  multivaluedness  and  approaching  zero  as  C  ♦  0 
more  slowly  than  any  definite  function.  Similarly,  the  coefficient  functions 
p(z)  in  (1)  here  a<initted  can  be  of  great  complexity,  especially  when  several 
irregular  singular  points  are  present,  and  a  useful  global  description  appears 
unlikely  in  the  general  case.  locally,  however,  the  class  of  potential 
functions  p(z)  can  be  described  by 

z'1  J^[p(t)/p(z)]1/2  dt  *  1  -  2Y  as  z  ♦  0  . 

As  indicated  in  the  Introduction,  the  conceptual  key  to  connection  lies 
in  the  two-variable  structure  of  the  Schroedinger  equation  emerging  from  (2) 
to  (4).  It  should  not  have  surprised  us  as  much  as  it  did,  for  it  is  already 
apparent  in  (1)  that  the  independent  variable  plays  two  distinct  physical 
roles.  The  first  term  in  (1)  represents  the  oscillatory  mechanism  and  its 
independent  variable  is  clearly  z/e,  with  local  wavelength  as  natural  unit, 
prompting  the  transformation  to 
(7)  x-i/*/p(t)dt  . 

By  contrast,  p(z)  represents  the  potential  and  its  variation,  which  is  not 
dependent  on  the  presence  of  waves,  and  the  role  of  z  in  it  is  therefore  a 
different  one.  The  formulation  sketched  in  this  Section  adds  the  insight  that 
the  dependence  on  the  modulation  variable  z  ♦♦  £x  »  5  enters  into  the 
Schroedinger  equation  (3)  only  through  a  relatively  minor  term 

♦(Ex)  =  xr'/r  -  y 

in  the  modulation  function  r'/r.  This  function  ♦  has  been  called 
"irregularity  function"  in  (IPM]  because  ♦  *  0  characterizes  the  regular 
singular  points. 
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The  critical  role  of  this  two-variable  structure  will  emerge  in  Section  4 


which  is  devoted  to  a  proof  that  connection  across  singular  points  is  a 
mathematical  process  local  in  ex,  even  though  asymptotic  in  x.  This  is, 
perhaps,  the  main  new  insight  gained  by  extending  the  fundamental  view  of 
connection  of  Olver  [1974]  to  irregular  singular  points.  It  explains  why  a 
merely  local  definition  of  $(ex)  -  -  and  thereby,  of  the  potential  and  of  the 
Schroedinger  equation  -  -  on  a  Riemann  surface  element 

A  =  {ex:  <  arg(ex)  <  2k,  0  <  |ex|  <  E  for  some  E  >  0} 

turns  out  sufficient.  (For  notational  convenience,  E  is  adjusted  so  that 
♦(ex)  is  analytic  up  to  the  rim  of  the  element  and  hence,  bounded  on  A.) 

In  a  shortwave  limit  e  -*■  0,  that  entails  little  restriction  on  the 
corresponding  domain  D  of  x.  On  account  of  the  two-variable  structure, 
moreover,  a  shortwave  limit  must  be  a  limit  ex  ♦  0.  Hindsight,  of  course, 
makes  all  of  this  appear  foreshadowed  in  the  structure  of  (1),  where  the 
first,  oscillation  term  is  defined  globally  in  z/e,  even  if  the  potential 
p(z)  be  defined  only  locally. 


3 .  Branch  Structure 

It  will  help  to  summarize  now  the  result  of  [IPM]  used  in  the  proof  of 
connection  in  Section  4  and  to  indicate  their  motivation.  If  a  limit  e  ♦  0 
is  taken  at  fixed  x,  then  t(ex)  +  0,  by  (5),  and  (3)  approaches  a  form  of 
Bessel's  equation  (which  observation  started  Langer's  [1931]  work).  Its 

singular  point  is  regular  with  Frobenius  exponents  0  and  1-2Y  >0.  If 

1  -2Y 

1-2Y  is  not  an  integer,  that  implies  solutions  f_(x)  and  x  f  (x)  with 

53  m 

entire  functions  fg,  fm  (and  ( 0 )  =  fm(0)  =  1),  which  turn  out  to  depend 

2  1 

only  on  x  .  Integer  values  of  —  -  Y  correspond  to  the  Frobenius  exceptions 

for  which  only  f„(x)  is  entire. 

J  m 
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It  is  plausible  that  (3)  may  have  an  analogous  fundamental  system 

w^en  e  ^  0#  which  displays  the  branch  structure  of  the  irregular 

singular  point  most  clearly.  If  ym/ys  *  0  as  x  -  0,  then  y8  has  there 

the  stronger  singularity  and  it  appears  natural  to  call  it  the  stronger 

solution  and  ym,  the  milder .  Such  a  system  has  been  constructed  [IPM]  to 

obtain  a  representation  of  the  branch  structure  at  the  general  irregular  point 

of  wave  modulation  and  in  particular,  to  find  out  what  replaces  the  entire 

functions  fs,  fm  and  to  explore  how  departure  from  entirety  can  be 

characterized.  The  underlying  idea  emerges  most  simply  in  the  following 

construction  of  yQ(x)  for  ^  >  Re  T  >  - 

s  2  2 

Since  (3)  can  be  written  (r^y1 )•  *  r^y,  a  simple  Volterra  equation 
associated  with  it  is 


(8) 


y*(x)  =  [r(x)J~2  /*fr(v)]  2y(v)dv 
y(x)  =  1  +  /*  y’ (v)dv  . 


By  (5),  a  simple  iterative  approach  to  (8)  is  by  a  sequence  (b  (x)}  such 

n 

that 


fx  p(ev)  2  fV'i  K  , 

0  [FT^O]  lxJ  bn(v,dv  ' 


bo  E 


b  (0) 
n 


for  n  >  1 


and  since  P(S>  is  a  mild  function,  it  can  be  estimated  at  the  expense  of  a 

small  power.  It  emerges  readily  in  this  way,  by  estimation  of  b*  from  (9) 

n 

and  (6)  and  thence,  bn ,  recursively,  that 

b  (x)  =  B  (ex) (x/2)2"  , 

n  n 

MO)  |0n(?)l  *  kn  “  . 

s  -  -Re  Y  -  j  ♦  lub  |$(u£) |  . 

u@  ( 0 , 1 ) 


Therefore,  if  |£x|  <  E(Y)  chosen  to  assure  s  <  0,  then  the  rapid  decrease 

of  k^  with  n  documents  a  majorant  series  assuring  the  convergence  of 

1  b  to  a  solution  y  (x)  of  (3)  analytic  on  the  Riemann  sector  D. 

Since  yo(0)  »  1,  J  b  (x)  generalizes  Frobenius'  entire  function  f„(x),  but 
s  n  3 

the  'coefficients'  8^  are  generally  multi-valued  functions  of  ex. 

Observe  that  (10)  suggests  y  **  J  b  tends  to  an  even  function  of  x 

s  n 

in  some  sense  as  ex  ♦  0.  This  can  be  made  more  precise  by  applying  the  same 

approach  to  the  estimation  of  |b'(x)  +  b'(x  exp-*i)|  and  thence, 

n  n 

|b  (x)  -  b  (x  exp-^i) |  to  show  [IPM]  that 
n  n 

|b  (X)  -  b  (xe‘Wl)|  <  6  (|ex|)n  k * | x/2 | 2" 
n  n  s  n 

where 

(ID  «  (|5|)  +  0  as  |C|  +  0  . 

s 

These  bounds  still  decrease  fast  enough  with  n  to  be  summed  to  a  bound  on 
the  degree  of  oddness  of  yg(x)  in  terms  of  the  modified  Bessel  function 
Iy(z)  [Olver  1974,  p.  60]:  For  >  Re  Y  >  -  j  and  x,  x  exp(-iri)  in  D 
and  | e x |  <  E(Y), 

(12)  |y  (x)  -  y  (xe_1,1)|  <  6  (  |  ex  | )  r(-s)  |  x/2 1  2+®I  (|x|)  . 

s  s  s  — s 

Thus  yg(x)  approaches  evenness  as  |ex|  ♦  0  uniformly  on  compact  subsets  of 

the  cut  x-plane.  For  fixed  ex,  on  the  other  hand,  yg(x)  may  lose  its 

evenness  exponentially  fast  as  |x|  increases. 

A  good  representation  of  the  milder  solution  ym(x)  of  (3)  depends  on 

1-2Y 

identification  of  the  exact  function  generalizing  the  factor  x  of 

fm(x)  .  It  turns  out  to  be  just  the  function  z(x)  defined  by  (2),  indeed 
(IPM] , 

(13)  z(x)  =  x1-2lr  C (ex) 
with  a  mild  function  C(?)  such  that 

(14)  ( 1 -2Y )p2?  ♦  -ie  as  5  ♦  0  . 
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Then  ym(x)/z(x)  =  y(x)  satisfies  a  differential  equation  related  to  (3)  and 
can  be  constructed  for  all  Re  Y  <  by  an  iteration  paralleling  that  just 
sketched  to  obtain  [IPM]  a  representation 

00 

(15)  y  (x)  =  z(x)[1  +  l  a  (Ex)(x/2)2n] 

ro  1  n 

with  bounds 

|a  (?)|  <  k  =  r (m)/[nl r (m+n) ]  , 

n  n 

where  m  =  3/2  -  Re  Y  -  6  (|ex|)  >  0  for  |£x|  <  another  E(Y),  since 

<5  (  | C I )  ♦  0  as  |S|  ♦  0.  Of  course,  C<5)  and  a  (£)  are  generally 
/.  n 

multivalued  functions,  but  «  (0)  is  defined  and  nonzero,  so  that 

n 

ym/z  »  y(x)  also  tends  to  an  even  function  as  Ex  ♦  0.  An  oddness  bound  is 

obtained  by  an  estimate  paralleling  that  indicated  above: 

Theorem  4  [IPM]  .  For  x  and  xe”"*  in  D  and  |ex|  <  E(Y) , 

ly(x)  -y(xe-1Ti)|  <  6  ( |Ex|  >r(m)  |x/2|2-mi  (|x|) 
in  m 

and  S  ( | 5 | )  ♦  o  as  |5|*0. 

m 

The  same  comment  therefore  applies  to  ym/z  »  y  as  follows  (12)  . 

No  similarly  simple  approach  has  been  found  yet  for  the  stronger 
solution  yg  for  Re  Y  <  -  where  the  simple  Volterra  equation  (8)  can ,  by 

(5),  be  used  only  at  the  price  of  a  regularization  of  the  first  integral  in 
(8) .  A  stronger  solution  in  the  sense  indicated  is  defined  only  up  to  an 
additive  multiple  of  the  milder  solution,  which  is  undesirable  for  a 
fundamental  system  displaying  clearly  the  branch  structure  of  the  singular 
point.  The  regularization  adopted  in  flPM]  avoids  that  additive  multiple  and 
constructs  a  stronger  solution  of  the  form 

00 

(16)  y  ( x)  -  1  +  I  B  (ex)(x/2)2n 

s  i  n 


-12- 


for  all  Re  Y  <  4,  but  the  estimates  of  |S  |  are  more  laborious,  and  less 
2  n 

close,  than  for  4  >  Re  Y  >  -  — .  Indeed,  8  (0)  is  found  to  exist  for  all  n 
2  2  n 

only  when  -  Re  Y  is  not  an  integer  and  the  near-evenness  of  yg(x)  holds 
only  for  them: 

Theorem  5  [IPM] .  For  non-integer  ^  -  Re  Y  >  0,  x  and  xe 

in  D  and  sufficiently  restricted  |£x|, 

ly_(x)  -  y  (xe"Wi)|  <  C(Y)«  ( |ex | > |x/2| 2+SI  (|x|), 

3  S  S  “S 

with  s  =  -ReY-4+  lub  |  <t>  <uC )  |  >  0  and  6  ( 1 5 1 )  ♦  0 

ue(0,1)  3 

as  |C|  ♦  0. 

4 .  Connection 

If  y(x)  satisfies  (3),  then  W(x)  =  r(x)y(x)  satisfies 

(17)  W  -  (1  +  r"/r)W 
and  by  ( 4 ) 

(18)  r"/r  =  x“2[Y(Y-1)  +  $(2Y-1  +  4»  + 

so  that  |r"/r|  is  integrable  along  paths  in  D  bounded  away  from  x  =  0. 
This  confirms  [Olver  1974,  p.  222]  existence  of  a  fundamental  "WKB”  solution 
pair 

(19)  W+(x)  -  a ( x ) ex ,  W_(x)  =  b(x)e_x 

with  functions  a(x),  b(x)  analytic  on  D  and  bounded  for  large  |x| 
(provided,  of  course,  e*  *  5  e  4  so  that  $  and  are  bounded).  This 

is  the  fundamental  system  of  (1)  most  strikingly  describing  the  asymptotic 
wave  character  (undamped  on  the  lines  where  x  is  pure  imaginary)  of  the 
solutions . 

The  "amplitude  functions"  a,  b  are  determined  only  up  to  a  constant 
factor,  but  apart  from  that,  the  decay  of  |r"/r|  at  large  |x|  suffices 
(Olver  1974,  p.  223,  224]  to  assure  limits  for  a(x)  and  b(x)  as  |x|  ♦  • 


with  atg  x  any  integer  multiple  of  if.  Those  limits  are  the  well-known  wave 


amplitudes  in  the  first-order  WKB-approximations  to  the  solutions  of  (1). 

Since  any  solution  of  (17)  must  be  a  linear  combination  of  W+  and  W_, 

(20)  r(x)y  (x)  =  a  (x)eX  +  b  (x)e  * 

in  in  in 

and  the  same  holds  with  subscript  s  in  place  of  m,  and  supposing  they  can 

■N.  %i 

be  normalized  satisfactorily,  then  a  ,...,b  jure  similarly  analytic  and 

m  s 

bounded  for  large  |x|  .  Since  the  lefthand  side  of  (20)  has  been  seen  in 

Section  3  to  be  multivalued,  not  all  of  a  ,...  b  can  be  entire,  and 

s  m 

symmetry  makes  it  plausible  that  all  of  them  will  usually  turn  out  to  be 

multivalued.  This  prompts  the  question 

-  '  2fti 

a  (•)  -  a  («e  )  =  ? 

m  in 

which  is,  in  fact,  a  connection  question  for  WKB  coefficients  [Olver  1974, 
p.  481]  . 

In  view  of  the  many  contexts  in  which  connection  is  important,  it  is 

natural  that  many  different  forms  of  the  connection  problem  are  found  in  the 

literature,  but  most  of  them  can  be  related  to  each  other  with  little  work, 

and  a  treatise  on  connection  for  simple  turning  points  is  found  in  Chapter  13 

of  [Olver  1974].  In  any  case,  the  problem  turns  on  relating  the  respective 

limits  which  represent  the  WKB  coefficients  on  different  domains,  and  when  it 

is  recognized  that  those  domains  correspond,  in  the  frame  of  the  natural 

variable,  to  sheets  of  the  Riemann  surface  of  the  solution,  the  form  of  the 

connection  question  just  arrived  at  is  seen  to  be  a  natural  one. 

By  contrast  to  the  functions  a(x),  b(x)  first  mentioned,  a  (x)  and 

in 

>■  * 
b^(x)  are  normalized  implicitly  by  the  normalization  of  ym(x)/z(x)  »  y(x) 

and  this  turns  out  to  introduce  an  e-dependence  into  the  normalization  of 

a  ,  b  .  For  fixed  £  i<  0,  moreover,  |x|  is  bounded  by  E/e  on  the 
m  re 

Riemann  sector  D  on  which  the  differential  equation  (1)  has  been  defined. 
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and  the  connection  question  can  therefore  be  posed  only  in  the  limit  £  ♦  0 


This  aspect  is  discussed  in  the  Appendix,  where  it  is  shown  that  the  functions 

a  /(PC)  =  a  (x;£x)  and  b  /(PC)  =  b  (x;£x)  , 

mm  mm 

N.  -V 

rather  than  a  and  b  themselves,  are  certain  to  have  limits  as  £  +  0 
m  m 

and  |x|  ♦  ®.  For  an  assuredly  meaningful  connection  question,  we  should 
therefore  rewrite  the  identity  (20)  as 

(21)  r(x)[P(5)C(C)]_1y  (x)  =  a  (x)  eX  +  b  (x)  e'X 

mm  m 

(with  explicit  mention  of  the  dependence  of  a,,,,  on  £  =  £x  omitted  to 

2*i 

focus  attention  now  on  arg  x)  and  ask  a  (®e  )  -  a  (“)  =  ? 

mm 

Since  (21)  is  an  identity  in  x  on  D,  it  holds  equally  at 

x  exp(-xi),  if  that  point  is  also  in  D.  If  exp(-xi)  be  abbreviated  by 

-  1-Y 

j,  then  since  ym  =  zy  and  rz  =  x  PC,  by  (5)  and  (13),  the  identity 
[y(x)  -  y(jx)]x1-Ye-|x|  =  (am(x)  -  jY_1bm( jx)lex_,x| 

“  tbm(x)  -  jY-1am(jx)]e'X"|x| 

also  holds  on  D.  Now  let  |£x|  ♦  0,  but  |x|  *  m  in  such  a  way  that  the 
lefthand  side  of  (22)  still  tends  to  zero.  That  this  does  indeed  define  a 
non-empty  set  of  "intermediate  limits",  in  the  terminology  of  singular- 
perturbation  theory  [Eckhaus  1979],  is  a  corollary  of  Theorem  4  because  lOlver 
1974,  p.  435] 

I  ( |  x  | )  '  |2Wx|"1/2e*X'  as  |x|  +  * 
m 

and,  e.g.,  |x|  »  | log  6  (|£x|)|  will  serve. 

m 

For  the  choices  arg  x  =  w  and  arg  x  *  2f,  respectively,  such  a  limit 


of  (22)  yields 

(23) 
whence 

(24) 


Fi  Y-1  2Fi  Y-1  xi 

b  (®e  )  -  j  a  («),  a  (“e  )  «  j  b  (®e  ) 

m  mm  m 


am( 


2»i 

"e  )  -  a  (- 


.  *  i« 


)  -  2i  sin(YF)  b  («e  )  , 

m  m 
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The  choice  arg  x  =■  0  adds 


y-1  -Hi 

(25)  a  (*)  -  jT  b(*e  ) 

re  m 

to  (23),  whence  the  answer  to  the  connection  question  for  is 

(26)  b  (“e"1)  -  b  (“e”*1)  =  2i  sin(Y")  a  (•)  . 

rare  in 

For  non-integer  j  -  Re  Y,  a  parallel  argument  for  the  stronger  solution 

starts  from  the  identity 

r(x)y  (x)  «  a  (x)eX  +  b  (x)e  X 
s  s  s 

analogous  to  (20)  to  deduce  that  the  normalization  yg(0)  “1#  yg(0)  “  0 
assures  boundedness  of 

a  (x)/P (£)  =*  a  (x)  and  b  (x)/p(£)  -  b  (x) 

8  8  S  S 

as  £  ♦  0  and  leads  via  the  identity  (22)  with  y,  m  and  1-Y  replaced, 

respectively,  by  ys,  s  and  Y ,  by  the  help  of  Theorem  5  to  the  same 

connection  formulae  (24),  (26)  for  a_,  b_  in  the  place  of  a  .  b  ,  because 

sin[(1-Y)*]  »  sin(Y»).  [It  is  this  independence  of  subscript  which  makes 

(24),  (26)  more  convenient  for  present  purposes  than  various  other  relations, 

such  as  a  ("  exp  2*i)  »  a  (•)  exp(-2Y*i),  also  implied  by  (23)  and  (25),  or 
n  m 

their  counterparts  for  as,  bg  obtained  by  replacement  of  1-Y  by  Y .] 

With  appropriate  interpretation,  moreover,  the  same  connection  formulae 

'  '  '  2-1 
relate  a  ,  b  and  a  ,  b  ,  respectively,  because  (p  C)  tends  to  a 
vn  m  8  8 

definite  limit  as  £  ♦  0,  by  (14),  and  by  (6), 


Since  any  solution  y(x)  of  (3)  is  a  linear  combination  of  the  fundamental 
pair  (yg,  ym) ,  the  functions  a(x)  and  b(x)  in  the  representation 

r(x)y(x)  -  a(x)ex  +  b(x)e~x 

in  terms  of  the  fundamental  pair  (W+,  W_)  of  (17)  are  linear  combinations  of 
a  ,  a  and  b  ,  b  ,  respectively,  and  therefore  satisfy  (24)  and  (26)  as 

a  III  8  re 

well.  In  the  limit  e  ♦  0  and  with  interpretation  appropriate  to  the 
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normalization  of  y(x),  the  connection  formulae  (24)  and  (26)  are  therefore  a 
general  corollary  of  (3)  under  the  two  hypotheses  of  Section  2,  at  least,  as 
long  as  Re  Y  -  j  is  not  a  negative  integer. 


Appendix 


The  somewhat  delicate  issue  of  normalization  for  connection  may  be 
brought  under  control  in  two  steps.  For  fixed  e,  the  domain  D  of  x  on 
which  r(x),  and  hence  also  the  differential  equation  (17),  is  defined  is  a 
Riemann  sector  of  radius  E(Y)/e.  Let  the  particular  functions  a,  b  in  ( 1 9) 
normalized  in  the  manner  of  Olver  [1974.  •  -s..  220-222]  be  denoted  by  a0, 

bQ .  Then  ag[bg]  =  1  and  ]  *  0  arf;  *  ->  lot  of  minD[raaxD]Re  x,  which 

depends  on  e,  and  thus  aQ  =  a g(x;£),  i:.  -  ,;,^(xf,e),  and  the  first  step  will 
be  to  confirm  that  their  dependence  cn  t  ..'akens  as  e  ♦  0 . 

Since  <MC)  is  analytic  and  bounded  on  the  Riemann  sector  A  for 
0  <  | Cl  <  E(Y),  the  same  follows  for 

x2*(x)  =  Y  (Y-1 )  +  (2Y-1  )♦(?)  +  *2  +  C<t>'  (C) 

=  x2r"/r 

in  (18)  because  it  tends  to  Y(Y-1)  as  C  *  ex  +  0,  by  (6) .  For  the  basic 
connection  question  of  Section  4,  it  is  sufficient  to  restrict  the  Riemann 
sector  D  to  a  disc  of  the  same  radius  cut  along  the  positive  real  axis  of 
x  for  a^(x;£),  and  along  the  negative  real  axis,  for  bg(x»e).  Olver's 
[1974,  p.  221]  variation  function  for  ag  is 

'|,(X)  "  1 E/e | 

evaluated  along  progressive  paths,  and  if  such  a  path  keeps  distance  R  from 

2 

the  origin,  then  since  x  is  bounded, 

(A1 )  ¥  -  0(R_1)  as  R  ♦  » 

The  same  holds  for  the  variation  function  for  bg,  which  differs  only  in  that 
the  lower  limit  is  +  |E/e|.  The  functions  furnish  [Olver  1974,  p.  221] 
bounds 

( A2 )  |a0(x>e)  -  1 | ,  |a^(xje)|  <  ey(x)  -  1  -  0(R_1) 

and  similarly,  for  bQ .  If  now  |e^|  <  |  e  |,  then  D(e^)  3  D(e  )  and  on 
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D(e^)  the  solution  W+  normalized  for  D(e^)  is  a  linear  combination  of 
W +,  W_  normalized  for  Dfe^),  thus 

a0(x;e.)eX  =  a0(x;ek)eX  +  Vx'ek>e’X 
with  constants  c^k,  d^  computed  from  the  respective  normalization  and 
bounds  to  yield 

-2(x+E/£  ) 

a0(x;Ei)  =  aQ(x;ek){l  +  O^)}  +  b^x.-e^e  Ofe^ 

and  on  DfED,  Re(x  +  E/£D  >0.  As  long  as  |x|  is  well  bounded  away  from 
O,  therefore,  agfxjE^  "  a0^X'Ek^  +  O(e^)  as  lEil  +  °#  and  similarly  for 
b0(x;e),  and  by  the  bounds  (A1),  (A2),  ag  and  bg  tend  to  limits  as 
I x  |  ♦  “  in  D (0  )  . 

For  the  second  step,  note  that  the  amplitude  functions  a  and  b  in 

m  m 

(20)  are  renormalizations  of  ag  and  bg  so  that 

a  (x;E)  =  A  a  (x;£),  b  < x;£)  =  b  b  (x»e) 
m  0  ra  o 

with  coefficients  A,  B  possibly  dependent  on  e.  From  !20),  (5)  and  (13), 
therefore, 

(A3)  x1  Yym(x)/z(x)  =  x1  Yy(x)  *  (PC )  ""^  [Aa^x; E)eX  +  BbQ(x;e)e  X]  . 

The  differential  equation  for  y(x)  is,  by  (2)  and  (3), 

y"  +  2(r'/r  +  z'/z)y'  *  y 

and  since  the  normalization  to  y(0)  »  1,  y' (0)  *  0  recognized  in  (15)  is 
independent  of  £,  y  inherits  from  (r'/r  +  z'/z)x  the  property  that  it 
depends  on  e  only  through  £  =  Ex.  Like  (r'/r  +  z'/z)x,  moreover,  y 
depends  continuously  on  £  in  A,  for  fixed  x  f>  0,  and  as  £  *  0,  by  (4) 
and  (13),  r'/r  +  z'/z  ♦  (1-Y)/x  and  the  differential  equation  for  y 
becomes  a  form  of  Bessel's,  with  solution 

( A4)  lim  y(x)  =  T(|  -  y)  (x/2 ) Y"  Ih  v(x)  . 

£*0  2  /2  'Y 

1-v„ 

This  shows  x  y(x)  to  tend  to  a  well-defined  function  of  x  on  D(0)  as 
£  ♦  0,  and  since  a0(x;e)  and  b0(x;e)  have  been  shown  to  tend  to  limit 
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functions  on  D(0)  as  e  +  o,  if  follows  from  (A3)  that  the  functions 


A /(PC)  and  B/ (PC )  must  tend  to  limits  as  e  -*■  0  and  C  *  Oj  of  course, 
these  limits  might  depend  on  the  direction  of  approach,  which  is  determined 
by  arg  x . 

In  sum. 


m 

PC 


A  ,  .  m  B 

PC  a0(x>  )f  PC  “  PC  b0(x’  * 


where  A/(PC),  B/(PC)  have  limits  as  e  ♦  0  and  5  ♦  0,  while  the  limits  as 
e  ♦  0  of  aQ,  bQ  are  defined  on  D(0)  and  tend  there  to  limits  as 
|x|  ♦ 
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